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Abstract 



Jet production in PbPb collisions at a nucleon-nucleon center-of-mass energy of 
2.76 TeV was studied with the CMS detector at the LHC, using a data sample cor- 
responding to an integrated luminosity of 6.7 }ib^^. Jets are reconstructed using the 
energy deposited in the CMS calorimeters and studied as a function of collision cen- 
trality. With increasing collision centrality, a striking imbalance in dijet transverse 
momentum is observed, consistent with jet quenching. The observed effect extends 
from the lower cut-off used in this study (jet pj = 120 GeV/c) up to the statistical limit 
of the available data sample (jet pj ~ 210 GeV/c). Correlations of charged particle 
tracks with jets indicate that the momentum imbalance is accompanied by a soften- 
ing of the fragmentation pattern of the second most energetic, away-side jet. The dijet 
momentum balance is recovered when integrating low transverse momentum parti- 
cles distributed over a wide angular range relative to the direction of the away-side 
jet. 
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1 Introduction 

High-energy collisions of heavy ions allow the fundamental theory of the strong interaction — 
Quantum Chromodynamics (QCD) — to be studied under extreme temperature and density 
conditions. A new form of matter formed at energy densities above ~1 GeV/ fm'' is pre- 
dicted in Lattice QCD calculations [Sj. This quark-gluon plasma (QGP) consists of an extended 
volume of deconfined and chirally-symmetric quarks and gluons. 

Heavy ion collisions at the Large Hadron Collider (LHC) are expected to produce matter at en- 
ergy densities exceeding any previously explored in experiments conducted at particle acceler- 
ators. One of the first experimental signatures suggested for QGP studies was the suppression 
of high-transverse-momentum (px) hadron yields resulting from energy loss suffered by hard- 
scattered partons passing through the medium [6J. This parton energy loss is often referred 
to as "jet quenching". The energy lost by a parton provides fundamental information on the 
thermodynamical and transport properties of the traversed medium, which is now believed to 
be strongly coupled as opposed to an ideal gas of quarks and gluons (recent reviews: [|Zn81). 
Results from nucleus-nucleus collisions at the Relativistic Heavy Ion Collider (RHIC) ||9]-[T2| 
have shown evidence for the quenching effect through the suppression of inclusive high-px 
hadron production and the modification of high-px dihadron angular correlations when com- 
pared to the corresponding results in much smaller systems, especially proton-proton colli- 
sions. Preliminary results for fully reconstructed jets at RHIC, measured in AuAu collisions 
at ^/s^ = 200 GeV ||T3tiT6| , also hint at broadened jet shapes due to medium-induced gluon 
radiation. 

Studying the modification of jets has long been proposed as a particularly useful tool for prob- 
ing the QGP properties [IT'.'TS'I. Of particular interest are the dominant "dijets", consisting of 
the most energetic ("leading") and second most energetic ("subleading") jets. At leading order 
(LO) and in the absence of parton energy loss, the two jets have equal px with respect to the 
beam axis and are emitted very close to back-to-back in azimuth (Aq^^ijet = | ^'jeti — 9]et2 1 ~ tt). 
However, medium-induced gluon emission in the final state can significantly alter the energy 
balance between the two highest-px jets and may give rise to large deviations from Ag^dijet ~ ^• 
Such medium effects in nuclear interactions are expected to be much larger than those due to 
higher-order gluon radiation, which is also present for jet events in pp collisions. The study 
of medium-induced modifications of dijet properties can therefore shed light on the transport 
properties of the QCD medium formed in heavy ion collisions. 

The dijet analysis presented in this paper was performed using the data collected in 2010 from 
PbPb collisions at a nucleon-nucleon center-of-mass energy of ^/s^ = 2.76 TeV at the Compact 
Muon Solenoid (CMS) detector. The CMS detector has a solid angle acceptance of nearly An and 
is designed to measure jets and energy flow, an ideal feature for studying heavy ion collisions. 
A total integrated (PbPb) luminosity of 8.7 ]dh^^ was collected, of which 6.7 ]dh^^ has been 
included in this analysis. Recently, related results on a smaller data sample (1.7 }ih^^) have 
been reported by ATLAS Ill9l . 

Jets were reconstructed based on their energy deposits in the CMS calorimeters. In general, 
the jet quenching effect on partons traversing the medium with different path lengths will lead 
to modifications in the observed dijet energy balance due to a combination of two effects: the 
radiated energy can fall outside the area used for the determination of the jet energy, and the 
energy can be shifted towards low momentum particles, which will not be detected in the 
calorimetric energy measurement. Such unbalanced events are easy to detect visually even at 
the level of event displays, and numerous examples were in fact seen during the first days of 
data taking (e.g. Fig.[l]l. 
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2 Experimental method 




Figure 1: Example of an unbalanced dijet in a PbPb collision event at = 2.76 TeV. Plot- 

ted is the summed transverse energy in the electromagnetic and hadron calorimeters vs. rj 
and (p, with the identified jets highlighted in red, and labeled with the corrected jet transverse 
momentum. 



The data provide information on the evolution of the dijet imbalance as a function of both 
collision centrality (i.e., the degree of overlap of the two colliding nuclei) and the energy of 
the leading jet. By correlating the dijets detected in the calorimeters with charged hadrons 
reconstructed in the high-resolution tracker system, the modification of the jet fragmentation 
pattern can be studied in detail, thus providing a deeper insight into the dynamics of the jet 
quenching phenomenon. 

The paper is organized as follows: the experimental setup, event triggering, selection and char- 
acterization, and jet reconstruction are described in Section|2] Section|3]presents the results and 
a discussion of systematic uncertainties, followed by a summary in Section |4] 



2 Experimental method 

The CMS detector is described in detail elsewhere ||20| . The calorimeters provide hermetic 
coverage over a large range of pseudorapidity, < 5.2, where rj = —In [ tan(0/2)] and 6 is 
the polar angle relative to the particle beam. In this study, jets are identified primarily using 
the energy deposited in the lead-tungstate crystal electromagnetic calorimeter (ECAL) and the 
brass/scintillator hadron calorimeter (HCAL) covering \t]\ < 3. In addition, a steel/ quartz- 
fiber Cherenkov calorimeter, called Hadron Forward (HE), covers the forward rapidities 3 < 
\rj\ < 5.2 and is used to determine the centrality of the PbPb collision. Calorimeter cells are 
grouped in projective towers of granularity in pseudorapidity and azimuthal angle given by 
Afj X Acp = 0.087 X 0.087 at central rapidities, having a coarser segmentation at forward rapidi- 
ties. The central calorimeters are embedded in a solenoid with 3.8 T central magnetic field. The 
event display shown in Eig.jl] illustrates the projective calorimeter tower granularity over the 
full pseudorapidity range. The CMS tracking system, located inside the calorimeter, consists 
of pixel and silicon-strip layers covering \t]\ < 2.5, and provides track reconstruction down to 
pj ~ 100 MeV/c, with a track momentum resolution of about 1% at pj = 100 GeV/c. A set 
of scintillator tiles, the Beam Scintillator Counters (BSC), are mounted on the inner side of the 
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HF calorimeters for triggering ar\d beam-halo rejection. CMS uses a right-handed coordinate 
system, with the origin located at the nominal collision point at the center of the detector, the 
:t-axis pointing towards the center of the LHC ring, the y-axis pointing up (perpendicular to 
the LHC plane), and the z-axis along the counterclockwise beam direction. The detailed Monte 
Carlo (MC) simulation of the CMS detector response is based on GEANT4 [21] . 

2.1 Data samples and triggers 

The expected cross section for hadronic inelastic PbPb collisions at a/s^ = 2.76 TeV is 7.65 b. 



corresponding to the chosen Glauber MC parameters described in Section 2.3 In addition, 
there is a sizable contribution from large impact parameter ultra-peripheral collisions (UPC) 
that lead to the electromagnetic breakup of one, or both, of the Pb nuclei II22| . particles per unit 
of pseudorapidity, depending on the impact parameter. 

For online event selection, CMS uses a two-level trigger system: Level-1 (LI) and High Level 
Trigger (HLT). The events for this analysis were selected using an inclusive single-jet trigger 
that required an LI jet with pj > 30 GeV/c and an HLT jet with pj > 50 GeV/c, where neither 
Pt value was corrected for the pj-dependent calorimeter energy response discussed in Sec- 



tion [2.4.3 The efficiency of the jet trigger is shown in Fig. |2|(a) for leading jets with \r]\ < 2 as a 
function of their corrected pj. The efficiency is defined as the ratio of the number of triggered 
events over the number of minimum bias events (described below). The trigger becomes fully 
efficient for collisions with a leading jet with corrected pj greater than 100 GeV/c. 




Figure 2: (a) Efficiency curve for the HLT 50 GeV/c single-jet trigger, as a function of the cor- 
rected leading jet transverse momentum. Error bars shown are statistical, (b) Correlation be- 
tween the number of pixel hits and HF total energy for a single run containing 60 k minimum 
bias events, after selections as described in the text. 



In addition to the jet data sample, a minimum bias event sample was collected using coinci- 
dences between the trigger signals from the +z and — z sides of either the BSC or the HF. This 
trigger has an efficiency of more than 97% for hadronic inelastic collisions. In order to sup- 
press non-collision related noise, cosmic rays, double-firing triggers, and beam backgrounds, 
the minimum bias and jet triggers used in this analysis were required to fire in time with the 
presence of both colliding ion bunches in the interaction region. It was checked that the events 
selected by the jet trigger described above also satisfy all triggers and selections imposed for 
minimum bias events. The total hadronic collision rate varied between 1 and 210 Hz, depend- 
ing on the number of colliding bunches (between 1x1 and 129 x 129) and on the bunch intensity. 
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2.2 Event selection 

In order to select a pure sample of inelastic hadronic collisions for analysis, a number of offline 
selections were applied to the triggered event sample, removing contaminations from UPC 
events and non-collision beam backgrounds (e.g. beam-gas). Table [T] shows the number of 
events remaining after the various selection criteria are applied. First, beam halo events were 
vetoed based on the timing of the +z and — z BSC signals. Then, to veto UPC and beam-gas 
events, an offline HP coincidence of at least three towers on each side of the interaction point 
was required, with a total deposited energy of at least 3 GeV. Next, a reconstructed vertex was 
required with at least two tracks of px > 75 MeV/c, consistent with the transverse beam spot 
position and the expected collision region along the z-axis. Pinally, to further reject beam-gas 
and beam-scraping events, the length of pixel clusters along the beam direction were required 
to be compatible with particles originating from the primary vertex. This last selection is identi- 
cal to the one used for the study of charged hadron pseudorapidity density and px spectrum in 
7 TeV pp collisions Il23l . Pigure|2](b) shows the correlation between the total energy deposited 
in the HP calorimeters and the number of hits in the first layer of the silicon pixel barrel detec- 
tor after these event selections. A tight correlation between the two detectors is observed, with 
very few of the events showing HP energy deposits that deviate significantly (at any given 
number of pixel hits) from the expectations for hadronic PbPb collisions. This correlation is 
important to verify the selection of a pure collision event sample, and also to validate the HP 
energy sum as a measure of event centrality (Section |2.3| . 

Starting from inelastic hadron collisions based on the selections described above, the basic off- 
line selection of events for the analysis is the presence of a leading calorimeter jet in the pseudo- 
rapidity range of | // 1 < 2 with a corrected jet px > 120 GeV/ c. By selecting these leading jets we 
avoid possible biases due to inefficiencies close to the trigger threshold. Furthermore, the se- 
lection of a rather large leading jet momentum expands the range of jet momentum imbalances 
that can be observed between the leading and subleading jets, as the subleading jets need a 
minimum momentum of px > 35-50 GeV/c to be reliably detected above the high-multiplicity 
underlying event in PbPb collisions (Section [2. 4| . In order to ensure high quality dijet selection, 
kinematic selection cuts were applied. The azimuthal angle between the leading and sublead- 
ing jet was required to be at least 2n/3. Also, we require a minimum px of px,i > 120 GeV/c 
for leading jets and of px,2 > 50 GeV/c for subleading jets. No explicit requirement is made 
either on the presence or absence of a third jet in the event. Prior to jet finding on the selected 
events, a small contamination of noise events from uncharacteristic PCAL and HCAL detector 
responses was removed using signal timing, energy distribution, and pulse-shape information 
Il24ll25l . As a result, about 2.4% of the events were removed from the sample. 

2.3 Centrality determination 

Por the analysis of PbPb events, it is important to know the "centrality" of the collision, i.e., 
whether the overlap of the two colliding nuclei is large or small. In this analysis, the observable 
used to determine centrality is the total energy from both HP calorimeters. The distribution of 
the HP signal used in this analysis is shown in Pig. |3](a). The shape of the energy distribution 
is characteristic of all observables related to (soft) particle production in heavy ion collisions. 
The more frequent peripheral events with large impact parameter produce very few particles, 
while the central ones with small impact parameter produce many more particles because of 
the increased number of nucleon-nucleon interactions. 

The distribution of this total energy was used to divide the event sample into 40 centrality 
bins, each representing 2.5% of the total nucleus-nucleus interaction cross section. Because 
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Table 1: Event selection criteria used for this analysis. The percentage of events remaining after 
each criterion, listed in the last column, are with respect to the previous criterion (the event 
selection criteria are applied in the indicated sequence). 



Criterion 


Events remaining 


% of events remaining 


Jet triggered events (p^™" > 50 GeV/c) 


149 k 


100.00 


No beam halo, based on the BSC 


148 k 


99.61 


HF offline coincidence 


111k 


74.98 


Reconstructed vertex 


110 k 


98.97 


Beam-gas removal 


110 k 


99.78 


ECAL cleaning 


107 k 


97.66 


HCAL cleaning 


107 k 


99.97 


> 2 jets with pj > 35 GeV/c and \r]\ <2 


71.9k 


67.07 


Leading jet p^i > 120 GeV/c 


4216 


5.86 


Subleading jet pj^2 > 50 GeV/c 


3684 


87.38 


A(pi2 of 2 jets > 2n/3 


3514 


95.39 



of inefficiencies in the minimum bias trigger and event selection, the measured multiplicity 
distribution does not represent the full interaction cross section. MC simulations were used 
to estimate the distribution in the regions where events are lost. Comparing the simulated 
distribution to the measured distribution, it is estimated that the minimum bias trigger and 
event selection efficiency is 97 ± 3%. 

For the jet analysis, these fine-grained bins were combined into five larger bins corresponding 
to the most central 10% of the events (i.e., smallest impact parameter), the next most central 
10% of the events (denoted 10-20%), and further bins corresponding to the 20-30%, 30-50%, 
and 50-100% selections of the total hadronic cross section. 

Simulations can be used to correlate centrality, as quantified using the fraction of the total 
interaction cross section, with more detailed properties of the collision. The two most com- 
monly used physical quantities are the total number of nucleons in the two lead (•^'^^Pb) nuclei 
which experienced at least one inelastic collision, denoted Npart/ and the total number of binary 
nucleon-nucleon collisions, Ncou- 

The centrality bins can be correlated to the impact parameter, b, and to average values and 
variances of Npart and N^ou using a calculation based on a Glauber model in which the nucle- 
ons are assumed to follow straight-line trajectories as the nuclei collide (for a review, see [26J). 
The bin-to-bin smearing of the results of these calculations due to the finite resolution and fluc- 
tuations in the HF energy measurement was obtained from fully simulated and reconstructed 
MC events generated with the ampt event generator [27] . Standard parameters of the Woods- 
Saxon function used to model the distribution of nucleons in the Pb nuclei were used [28] . The 
nucleon-nucleon inelastic cross section, which is used to determine how close the nucleon tra- 
jectories need to be in order for an interaction to occur, was taken to be 64 ± 5 mb, based on a fit 
of the existing data for total and elastic cross-sections in proton-proton and proton-antiproton 
collisions [29|. The uncertainties in the parameters involved in these calculations contribute to 
the systematic uncertainty in Npart and NcoU for a given bin. The other source of uncertainty in 
the centrality parameters comes from the determination of the event selection efficiency. 

Using the procedure outlined above, the mean and spread (RMS) values of the impact param- 
eter, Npart/ and N^on for the five bins used in this analysis, and their systematic uncertainties. 
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Sum HF Energy (TeV) Centrality Bin 

Figure 3: (a) Probability distribution of the total HF energy for minimum bias collisions (black 
open histogram). The five regions correspond to the centrality ranges used in this analysis. 
Also shown is the HF energy distribution for the subset of events passing the HLT jet trigger 
(red hatched histogram), (b) Distribution of the fraction of events in the 40 centrality bins for 
minimum bias (black open histogram) and HLT jet triggered (red hatched histogram) events. 
The centrality-bin labels run from 100% for the most peripheral to 0% for the most central 
events. 



were extracted and are listed in Table |2j The RMS values for the centrality parameters are 
due to their correlation with the percentage cross section and the width of the chosen centrality 
bins. 

It is important to note that the selection of rare processes, such as the production of high-px 
jets, leads to a strong bias in the centrality distribution of the underlying events towards more 
central collisions, for which Ncou is very large. This can be seen in Fig. |3] (a), where the HF 
energy distribution for events selected by the jet trigger is shown in comparison to that for 
minimum bias events. The bias can be seen more clearly in Fig.|3](b), where the distribution of 
minimum bias and jet-triggered events in the 40 centrality bins is shown. 

Table 2: Mean and RMS values for the distributions of impact parameter, b, number of partici- 
pating nucleons, Npart, and number of nucleon-nucleon collisions, NcoU/ for the centrality bins 
used in this analysis. The RMS values represent the spread of each quantity within the given 
bins due to the range of percentage cross section included. 



Centrality 


b mean (fm) 


b RMS (fm) 


Npaii mean 


Npart RMS 


NcoU mean 


NcoU RMS 


0-10% 


3.4 ± 0.1 


1.2 


355 ±3 


33 


1484 ± 120 


241 


10-20% 


6.0 ± 0.2 


0.8 


261 ±4 


30 


927 ± 82 


183 


20-30% 


7.8 ± 0.2 


0.6 


187 ±5 


23 


562 ± 53 


124 


30-50% 


9.9 ± 0.3 


0.8 


108 ±5 


27 


251 ± 28 


101 


50-100% 


13.6 ± 0.4 


1.6 


22 ± 2 


19 


30 ±5 


35 
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2.4 Jet reconstruction in PbPb collisions 

2.4.1 Jet algorithm 

The baseline jet reconstruction for heavy ion collisions in CMS is performed with an iterative 
cone algorithm modified to subtract the soft underlying event on an event-by-event basis ||30| . 
Each cone is selected with a radius AR = A(p^ + Arj^ = 0.5 around a seed of minimum trans- 
verse energy of 1 GeV. The underlying event subtraction algorithm is a variant of an iterative 
"noise /pedestal subtraction" technique [31J. Initially, the mean value, (Eceii)/ arid dispersion, 
(7"(Eceii)/ of the energies recorded in the calorimeter cells are calculated for all rings of cells that 
have at least 0.3 GeV transverse energy deposit at constant pseudorapidity. The algorithm sub- 
tracts (Eceii) + ^{^ceii) from each cell. If a cell energy is negative after subtraction, the value is 
set to zero. Subtracting the mean plus the dispersion, as opposed to simply the mean, compen- 
sates for the bias caused by the "zeroing" of negative-energy cells. Jets are then reconstructed, 
using a standard iterative cone algorithm ||32l l33l , from the remaining cells with non-zero en- 
ergy. In a second iteration, the pedestal function is recalculated using only calorimeter cells 
outside the area covered by reconstructed high-pj jets (px > 10 GeV/c). The threshold of 
10 GeV/c was chosen in studies optimizing the final extracted jet pj resolution. The cell en- 
ergies are updated with the new pedestal function (again subtracting mean plus dispersion) 
and the jets are reconstructed again, using the updated calorimeter cells. The performance of 
this algorithm is documented in Ref. II30I . Jet corrections for the calorimeter response have 
been applied, as determined in studies for pp collisions Ii34| . When applying the algorithm to 
PbPb data, the subtracted background energy for R = 0.5 jet cones ranges from 6-13 GeV for 
peripheral events (centrality bins 50-100%) to 90-130 GeV for central collisions (0-10%), before 
applying jet energy scale corrections. 

To perform a cross-check of the main results, the anti-Zcj algorithm [35 1 with a resolution pa- 
rameter of 0.5 was used to reconstruct jets, as was done for the pp reference measurements pre- 
sented here. The energy attributed to the underlying event was estimated and subtracted using 
the "average energy per jet area" procedure provided by the fastjet package ||36ll37l . In order 
to eliminate biases in the underlying event estimation, an ;y-strip of total width Arj = 1.6 cen- 
tered on the jet position was used, with the two highest energy jets in each event excluded Il38| . 
In addition, the anti-fej jets were reconstructed based on particle flow objects ||39l |40l instead 
of calorimeter-only information. A good agreement was found with the calorimeter-based, 
iterative cone algorithm results. 

2.4.2 Simulated data samples 

For the analysis of dijet properties in PbPb events, it is crucial to understand how the jet recon- 
struction is modified in the presence of the high multiplicity of particles produced in the PbPb 
underlying event. The jet-finding performance was studied using dijets in pp collisions simu- 
lated with the PYTHIA event generator (version 6.423, tune D6T) [41], modified for the isospin 
content of the colliding nuclei [42]. In order to enhance the number of Pythia dijets in the mo- 
mentum range studied, a minimum pj selection of 80 GeV/c was used. Lower pj selections, 
as discussed in [43J, were also investigated and found to agree with the pj = 80 GeV/ c results 
within uncertainties. The PYTHIA dijet events were processed with the full detector simulation 
and analysis chain. Additional samples were produced in which the pythia dijet events were 
embedded into a minimum bias selection of PbPb events at the raw data level. For this em- 
bedding procedure, both real PbPb data events (pythia+data), and PbPb events simulated 
with the HYDJET event generator [|42l (pythia+hydjet) were used. The hydjet parameters 
were tuned to reproduce the total particle multiplicities at all centralities and to approximate 
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the underlying event fluctuations seen in data. The hydjet events included the simulation 
of hard-scattering processes for which radiative parton energy loss was simulated, but coUi- 
sional energy loss was turned off [42 J. Both embedded samples were propagated through the 
standard reconstruction and analysis chain. 

The PYTHIA+DATA sample was used in several ways for studies of calorimeter jets. First, by 
matching the same PYTHIA dijet event reconstructed with and without the PbPb underlying 
event, the degradation of the jet pj and position resolution, the jet pj scale, and the jet-finding 
efficiency were determined as a function of collision centrality and jet pj (Section 2.4.3| |. In 



addition, pythia+data events were compared to non-embedded pythia for dijet observables 
such as azimuthal correlations and momentum balance distributions. Finally, to separate ef- 
fects due to the medium itself from effects simply due to reconstructing jets in the complicated 
environment of the underlying PbPb event, a direct comparison of results for PYTHIA+DATA 
and actual data events was made (Section |3Jj. 

The PYTHIA+HYDJET sample was used for studies of track momentum balance and track-jet 
correlations (Sections |3.2| and [33| , where access to the full MC particle level (truth) information 
for charged tracks is important for systematic studies. 

2.4.3 Jet-finding performance 

A detailed characterization of the CMS calorimeter jet-finding performance in pp collisions 
can be found in [44J. The dependence of the jet energy scale and of the jet energy resolution 
on centrality was determined using the pythia+data sample (Fig. |4j standard pp jet energy 
corrections are applied). In this study, reconstructed jets were matched to the closest generator- 
level jet in rj-(p within a cone of AR = 0.3. The residual jet energy scale dependence and the 
relative jet energy resolution are derived from the mean and standard deviation of the ap- 
proximately Gaussian distributions of the ratio of the reconstructed calorimeter jet transverse 
momentum p^^'"^*^' and the transverse momentum of jets reconstructed based on event gener- 
ator level final state particles p^'^"-''^*. For peripheral events in the 50-100% centrality selection, 
the jet energies are under-corrected by 5% after applying the standard pp jet energy corrections, 
and the jet energy resolution is found to be about 15% worse than in pp collisions. For the most 
central events, the large transverse energy per unit area of the underlying event leads to an 
over-correction of low-px jet energies by up to 10% and a degradation of the relative resolution 
by about 30% to (7(p^'''°J''Vpt''"^''*) = 0.16 at px = 100 GeV/c. The effect of the underlying 
event on the jet angular resolution was also studied. Integrated over jet pj > 50 GeV/c, the 
angular resolution in (p worsens from 0.03 for peripheral events (50-100%) to 0.04 for central 
events (0-10%), while the resolution in rj changes from 0.02 to 0.03 over the same centrality 
range. 

The jet reconstruction efficiency as a function of jet pj and centrality was extracted from the 
pythia+data sample as well, with the results shown in Fig. |5] For peripheral events, a jet- 
finding efficiency of 95% was found for a jet pj = 50 GeV/ c, while for central collisions the 
efficiency drops to 88% at the same pj. Jets with pj > 70 GeV/c are found with an efficiency 
greater than 97% for all collision centralities. No correction for the inefficiency near the thresh- 
old was applied in the subsequent analysis, as the effects of the reconstruction inefficiency are 
included in the pythia+data reference analysis. 

Finally, the rate of calorimeter jets reconstructed from fluctuations in the underlying event with- 
out the presence of a fragmenting pj parton, so called fake jets, for the jet selection used in this 
paper was determined using fully simulated 0-10% central hydjet events. Reconstructed jets 
in this sample are classified as fake jets if no matching generator-level jet of px > 20 GeV/c 
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Figure 4: The top row shows the mean of the ratio of reconstructed to generated jet momenta, 
{pj^^°^'^^/ Pj^'^^'^^), as a function of Pj^"^*^^, while the bottom row shows the relative resolution, 
i.e., the standard deviation of pj^^°^'^^/ p^*^"-^*^'. The standard pp jet energy corrections are in- 
cluded Filled circles are for the leading jets and open squares are for the subleading 
jets. The left, center, and right columns are for jets in PYTHIA+DATA events with centrality 50- 
100%, 20-30%, and 0-10%, respectively. On the jet resolution plots (bottom row), the dashed 
line is a fit to the leading jet resolution in pp events. The vertical bars denote the statistical 
uncertainty. 



is found within an r]-<p distance to the reconstructed jet axis smaller than 0.3. For leading jets 
with pt,i > 120 GeV/c, a fake jet fraction of less than 0.02% is found. In events with a pj^ > 
120 GeV/c leading jet, the fake jet fraction on the away-side of the leading jet {Acpu > 2tz/3) 
is determined to be 3.5% for reconstructed jets with pT,2 > 50 GeV/c and less than 0.02% for 
Pt^2 > 120 GeV/c. The effects of the degradation of jet performance in terms of energy scale. 



resolution, efficiency, and fake rate on the dijet observables are discussed in Section 3.1 



3 Results 

The goal of this analysis is to characterize possible modifications of dijet properties as a function 
of centrality in PbPb collisions. In addition to the standard event selection of inelastic hadronic 
collisions and the requirement of a leading jet with px^i > 120 GeV/c (Section |2.2| , most of the 
subsequent analysis required the subleading jet in the event to have pT,2 > 50 GeV/c, and the 
azimuthal angle between the leading and subleading jet (A(^i2) to be larger than 2n/3. Only 
jets within |?/| < 2 were considered for the analysis of calorimeter jets in Section 3.1 For a 
data set of Lint = 6.7 }ih^^, this selection yields 3514 jet pairs. For studies of correlations of 
calorimeter jets with charged particles (Sections |3.2| and [33] l, a more restrictive pseudorapidity 
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Figure 5: Jet reconstruction efficiency as a function of generator level jet pj for the leading 
jet (filled circles) and subleading jet (open squares). From left to right three centrality bins are 
shown: 30-100%, 10-30%, 0-10%. The vertical bars denote the statistical uncertainty. 



selection was applied. The analysis was performed mostly in five bins of collision centrality: 
0-10%, 10-20%, 20-30%, 30-50%, and 50-100%. 

Thus far, no pp reference data exist at the PbPb collision energy of ^/s^ = 2.76 TeV. Through- 
out the paper, the results obtained from PbPb data will be compared to references based on the 



PYTHIA and PYTHIA+DATA samples described in Section 2.4.2 



For most results, the pythia+data events will be used for direct comparisons. To calibrate 
the performance of pythia for the observables used in this analysis, the dijet analysis was also 
performed using the anti-A:T algorithm on 35 pb^^ of pp data at ^/s = 7 TeV, collected by CMS 
prior to the heavy ion data taking and compared to pythia simulations for the same collision 
system and energy. The same jet selection criteria used for the 2.76 TeV PbPb data were applied 
to both pp data and pythia. 

3.1 Dijet properties in pp and PbPb data 

The correlation between the transverse momentum of the reconstructed leading and sublead- 
ing jets in the calorimeters is plotted in Fig. |6] The top row contains PbPb data for peripheral, 
mid-central, and central events, the second row shows pp jets simulated by pythia and em- 
bedded into PbPb data, and the bottom panel shows pp jets from pythia without embedding. 
One can already observe a downward shift in the subleading jet pj for the more central PbPb 
events. In the following discussion, a more quantitative and detailed assessment of this phe- 
nomenon will be presented. 

3.1 .1 Leading jet spectra 

Figure [7| (a) shows the leading jet pj distributions for 7 TeV pp data and corresponding pythia 
simulations. The distribution of leading jet pj for PbPb is shown in Figs. [7| (b)-(f) for five dif- 
ferent centrality bins. The spectra obtained for PbPb data are shown as solid markers, whereas 
the hatched histograms show the leading jet spectrum reconstructed from pythia+data dijet 
events. All spectra have been normalized to unity. The detector-level leading jet spectra in 
PbPb data and the corresponding results for pythia+data samples show good quantitative 
agreement in all centrality bins over the pj range studied. 



3.1 Dijet properties in pp and PbPb data 
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Figure 6: Subleading jet px vs. leading jet px distributions. The top two rows show results for 
centrality 30-100% (left column), 10-30% (middle column) and 0-10% (right column), for PbPb 
data (top row) and reconstructed pythia jets embedded into PbPb data events (middle row). 
The panel in the bottom row shows the distribution for reconstructed jets from PYTHIA alone. 

It is important to note that the jet momentum spectra at detector level presented here have 
not been corrected for smearing due to detector resolution, fluctuations in/ out of the jet cone, 
or underlying event fluctuations. Therefore, a direct comparison of these spectra to analytical 
calculations or particle-level generator results is not possible. For the jet asymmetry and dijet 
A(p distributions discussed below, the effect of the finite jet energy resolution is estimated using 
the PYTHIA+DATA events. 

3.1.2 Dijet azimuthal correlations 

One possible medium effect on the dijet properties is a change of the back-to-back alignment of 
the two partons. This can be studied using the event-normalized differential dijet distribution, 
{l/N){dN/ dA(pi2), versus Acpu- Figure |8] shows distributions of Acpu between leading and sub- 
leading jets which pass the respective pj selections. In Fig. |8] (a), the dijet Acpu distributions 
are plotted for 7 TeV pp data in comparison to the corresponding pythia simulations using 
the anti-Zcx algorithm for jets based on calorimeter information, pythia provides a good de- 
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Figure 7: Leading jet pj distribution for dijet events with sub leading jets of Pt,2 > 50 GeV/c 
and A(pi2 > 2n/3 for 7 TeV pp collisions (a) and 2.76 TeV PbPb collisions in several centrality 
bins: (b) 50-100%, (c) 30-50%, (d) 20-30%, (e) 10-20% and (f) 0-10%. Data are shown as black 
points, while the histograms show (a) PYTHIA events and (b)-(f) PYTHIA events embedded into 
PbPb data. The error bars show the statistical uncertanties. 

scription of the experimental data, with slightly larger tails seen in the PYTHIA simulations. A 
recent study of azimuthal correlations in pp collisions at 7 TeV can be found in 1451 . For the 
PYTHIA comparison to PbPb results at ^/s^ = 2.76 TeV, this discrepancy seen in the higher 
energy pp comparison is included in the systematic uncertainty estimation. It is important to 
note that the PYTHIA simulations include events with more than two jets, which provide the 
main contribution to events with large momentum imbalance or A^i2 far from n. 

Figures|8](b)-(f) show the dijet Acpu distributions for PbPb data in five centrality bins, compared 
to PYTHIA+DATA simulations. The distributions for the four more peripheral bins are in good 
agreement with the pythia+data reference, especially for A^i2 ^ 2. The three centrality bins 
spanning 0-30% show an excess of events with azimuthally misaligned dijets {A(pi2 ^ 2), com- 
pared with more peripheral events. A similar trend is seen for the PYTHIA+DATA simulations, 
although the fraction of events with azimuthally misaligned dijets is smaller in the simulation. 
The centrality dependence of the azimuthal correlation in pythia+data can be understood as 
the result of the increasing fake-jet rate and the drop in jet reconstruction efficiency near the 
50 GeV/c threshold from 95% for peripheral events to 88% for the most central events. In PbPb 
data, this effect is magnified since low-px away-side jets can undergo a sufficiently large energy 
loss to fall below the 50 GeV/c selection criteria. 

Furthermore, a reduction of the fraction of back-to-back jets above A(/)i2 ^ 3 is observed for 
the most central bin. This modification of the A(pi2 distribution as a function of centrality can 
be quantified using the fraction Kg of dijets with A^i2 > 3.026 , as plotted in Fig. |9} for pT,i > 
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Figure 8: A(pi2 distributions for leading jets of pT,i > 120 GeV/c with subleading jets of pj^2 > 
50 GeV/c for 7 TeV pp collisions (a) and 2.76 TeV PbPb collisions in several centrality bins: (b) 
50-100%, (c) 30-50%, (d) 20-30%, (e) 10-20% and (f) 0-10%. Data are shown as black points, 
while the histograms show (a) PYTHIA events and (b)-(f) PYTHIA events embedded into PbPb 
data. The error bars show the statistical uncertainties. 

120 GeV/c and Pt,2 > 50 GeV/c. The threshold of 3.026 corresponds to the median of the 
A(^i2 distribution for PYTHIA (without embedding). The results for both the PbPb data and 
PYTHIA+DATA dijets are shown as a function of the reaction centrality, given by the number 



of participating nucleons, Npart/ as described in Section 2.3 This observable is not sensitive 
to the shape of the tail at A(pi2 < 2 seen in Fig. |8} but can be used to measure small changes 
in the back-to-back correlation between dijets. A decrease in the fraction of back-to-back jets 
in PbPb data is seen compared to the pure PYTHIA simulations. Part of the observed change 
in Rb (A^) with centrality is explained by the decrease in jet azimuthal angle resolution from 
dip = 0.03 in peripheral events to dtp = 0.04 in central events, due to the impact of fluctuations 
in the PbPb underlying event. This effect is demonstrated by the comparison of PYTHIA and 
PYTHIA+DATA results. The difference between the pp and pythia+data resolutions was used 
for the uncertainty estimate, giving the dominant contribution to the systematic uncertainties, 
shown as brackets in Fig. |9] 

3.1 .3 Dijet momentum balance 

To characterize the dijet momentum balance (or imbalance) quantitatively, we use the asym- 
metry ratio, 

Aj=P-^^, (1) 
Pt,i + Pt,2 
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Figure 9: Fraction of events with A(pi2 > 3.026 as a function of Npart/ among events with px,i > 
120 GeV/c and Pt,2 > 50 GeV/c. The result for reconstructed pythia dijet events (blue filled 
star) is plotted at Npart = 2. The other points (from left to right) correspond to centrality bins 
of 50-100%, 30-50%, 20-30%, 10-20%, and 0-10%. The red squares are for reconstruction of 
PYTHIA+DATA events and the filled circles are for the PbPb data, with statistical (vertical bars) 
and systematic (brackets) uncertainties. 



where the subscript 1 always refers to the leading jet, so that Aj is positive by construction. The 
use of Aj removes uncertainties due to possible constant shifts of the jet energy scale. It is im- 
portant to note that the subleading jet Pt,2 > 50 GeV/c selection imposes a px,! -dependent limit 
on the magnitude of Aj. For example, for the most frequent leading jets near the 120 GeV/c 
threshold, this limit is Aj < 0.41, while the largest possible Aj for the present dataset is 0.7 for 
300 GeV/c leading jets. Dijets in which the subleading jet is lost below the 50 GeV/c threshold 
are not included in the Aj calculation. 



In Fig. [10] (a), the Aj dijet asymmetry observable calculated by PYTHIA is compared to pp data 
at i/s = 7 TeV. Again, data and event generator are found to be in excellent agreement. This 
observation, as well as the good agreement between pythia+data and the most peripheral 



PbPb data shown in Fig. 10 (b), suggests that pythia at y/s = 2.76 TeV can serve as a good 



reference for the dijet imbalance analysis in PbPb collisions. 



The centrality dependence of Aj for PbPb collisions can be seen in Figs. 10 (b)-(f), in compar- 
ison to PYTHIA+DATA simulations. Whereas the dijet angular correlations show only a small 
dependence on collision centrality, the dijet momentum balance exhibits a dramatic change in 
shape for the most central collisions. In contrast, the PYTHIA simulations only exhibit a modest 
broadening, even when embedded in the highest multiplicity PbPb events. 

Central PbPb events show a significant deficit of events in which the momenta of leading and 
subleading jets are balanced and a significant excess of unbalanced pairs. The large excess of 
unbalanced compared to balanced dijets explains why this effect was apparent even when sim- 
ply scanning event displays (see Fig. [l]|. The striking momentum imbalance is also confirmed 
when studying high-px tracks associated with leading and subleading jets, as will be shown 
in Section [32| It is consistent with a degradation of the parton energy, or jet quenching, in the 
medium produced in central PbPb collisions. 
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Aj=(Pt,i-Pt/(Pt,i+Pt,2) 

Figure 10: Dijet asymmetry ratio, Aj, for leading jets of px,! > 120 GeV/c, subleading jets of 
Pt,2 >50 GeV/c and A(pi2 > 2n/3 for 7 TeV pp collisions (a) and 2.76 TeV PbPb collisions in 
several centrality bins: (b) 50-100%, (c) 30-50%, (d) 20-30%, (e) 10-20% and (f) 0-10%. Data are 
shown as black points, while the histograms show (a) PYTHIA events and (b)-(f) PYTHIA events 
embedded into PbPb data. The error bars show the statistical uncertainities. 



The evolution of the dijet momentum balance illustrated in Fig. 10 can be explored more quan- 
titatively by studying the fraction of balanced jets in the PbPb events. The balanced fraction, 
Rb(^/ < 0.15), is plotted as a function of collision centrality (again in terms of Npart) in Fig. 11 
It is defined as the fraction of all events with a leading jet having px4 > 120 GeV/c for which 
a subleading partner with Aj < 0.15 and Acpu > 2ti/3 is found. Since Rb{A] < 0.15) is cal- 
culated as the fraction of all events with px,! > 120 GeV/c, it takes into account the rate of 
apparent "mono-jet" events, where the subleading partner is removed by the px or A(p selec- 
tion. 

The Aj threshold of 0.15 corresponds to the median of the Aj distribution for pure PYTHIA 



dijet events passing the criteria used for Fig. 10 By definition, the fraction Rb{Aj < 0.15) of 



balanced jets in PYTHIA is therefore 50%, which is plotted as a dashed line in Fig. nil As will be 



discussed in Section 3.3 a third jet having a significant impact on the dijet imbalance is present 



in most of the large-Aj events in PYTHIA. 



The change in jet-finding performance from high to low px, discussed in Section 2.4.3 leads to 
only a small decrease in the fraction of balanced jets, of less than 5% for central PYTHIA+DATA 
dijets. In contrast, the PbPb data show a rapid decrease in the fraction of balanced jets with 
collision centrality. While the most peripheral selection shows a fraction of balanced jets of 
close to 45%, this fraction drops by close to a factor of two for the most central collisions. This 
again suggests that the passage of hard-scattered partons through the environment created in 
PbPb collisions has a significant impact on their fragmentation into final-state jets. 
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Figure 11: Fraction of all events with a leading jet with px,i > 120 GeV/cforwhichasubleading 
jet with Aj < 0.15 and A^i2 > In/S was found, as a function of Npart- The result for recon- 
structed PYTHIA dijet events (blue filled star) is plotted at Npart = 2. The other points (from 
left to right) correspond to centrality bins of 50-100%, 30-50%, 20-30%, 10-20%, and 0-10%. 
The red squares are for reconstruction of pythia+data events and the filled circles are for the 
PbPb data, with statistical (vertical bars) and systematic (brackets) uncertainties. 



The observed change in the fraction of balanced jets as a function of centrality, shown in Fig. 11 
is far bigger than the estimated systematic uncertainties, shown as brackets. The main contri 
butions to the systematic uncertainties include the uncertainties on jet energy scale and reso- 
lution, jet reconstruction efficiency, and the effects of underlying event subtraction. The uncer- 
tainty in the subtraction procedure is estimated based on the difference between pure PYTHIA 
and PYTHIA+DATA simulations. For central events, the subtraction procedure contributes the 
biggest uncertainty to Rb{Aj), of close to 8%. The uncertainty on the residual jet energy scale 
was estimated based on the results shown in the top row of Fig. |4j The full difference between 
the observed residual correction and unity, added in quadrature with the systematic uncer- 
tainty obtained for pp ||34| , was used as the systematic uncertainty on the jet pj and propagated 
to Rb{Aj). For the jet px resolution uncertainty, the full difference of the pythia+data result 
to the pp resolution, as shown in Fig. |4] (bottom), was used as an uncertainty estimate for the 
PbPb jet px resolution. The uncertainties in jet energy scale and jet resolution contribute 5% 
and 6%, respectively, to the 11% total systematic uncertainty in central events. For peripheral 
events, the total uncertainty drops to 9%, mostly due to the smaller uncertainty related to the 
PbPb background fluctuations for lower multiplicity events. 

3.1 .4 Leading jet pj dependence of dijet momentum imbalance 

The dependence of the jet modification on the leading jet momentum can be studied using the 
fractional imbalance Apx„,; = (Pt,i — Pt,2)/Pt,i- The mean value of this fraction is presented as 



a function of px,i in Fig. 12 for three bins of collision centrality, 30-100%, 10-30% and 0-10%. 
PYTHIA is shown as stars, pythia+data simulations are shown as squares, while the data are 
shown as circles. Statistical and systematic uncertainties are plotted as error bars and brackets, 
respectively. The dominant contribution to the systematic uncertainty comes from the observed 
px dependence of the residual jet energy correction in PbPb events (6% out of a total systematic 
uncertainty of 8%). The jet energy resolution and underlying event subtraction uncertainties 



3.2 Track-jet correlations 
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Figure 12: Mean value of the fractional imbalance (px,! — Pxi)/ Pt,i as a function of leading jet 
Pt for three centrality bins. The PbPb data are shown as circles with vertical bars and brack- 
ets indicating the statistical and systematic uncertainties, respectively. Results for PYTHIA are 
shown with blue stars, and pythia+data with red squares. The dot-dashed line to guide the 
eye is drawn at the value for pure PYTHIA for the lowest pj bin. 



The fractional imbalance exhibits several important features: the imbalance seen in PbPb data 
grows with collision centrality and reaches a much larger value than in PYTHIA or PYTHIA+DATA, 
In addition, the effect is clearly visible even for the highest-px jets observed in the data set, 
demonstrating that the observed dijet imbalance is not restricted to the threshold region in our 
leading jet selection. Within the present uncertainties, the px,i dependence of the excess imbal- 
ance above the PYTHIA prediction is compatible with either a constant difference or a constant 
fraction of px,i . 

The main contributions to the systematic uncertainty in (px,i — px,2)/pT,i are the uncertainties 
in the px-dependent residual energy scale (based on results shown in the top row of Fig. |4]|, 
and the centrality-dependent difference observed between PYTHIA and pythia+data seen in 
Fig. [12] As before, the uncertainty on the residual jet energy scale was estimated using the full 
difference between the observed residual correction and unity, and also assuming that within 
these limits the low-px and high-px response could vary independently. 

3.2 Track-jet correlations 

The studies of calorimeter jets show a strong change of the jet momentum balance as a func- 
tion of collision centrality. This implies a corresponding modification in the distribution of 
jet fragmentation products, with energy being either transported out of the cone area used to 
define the jets, or to low-momentum particles which are not measured in the calorimeter jets. 
The CMS calorimeter is less sensitive to these low momentum particles, or they do not reach 
the calorimeter surface. Information about changes to the effective fragmentation pattern as a 
function of Aj can be obtained from track-jet correlations. For this analysis, PYTHIA+HYDJET 
simulations are used as MC reference, to allow full access to MC truth (i.e., the output of the 
generator) information for tracks in the dijet signal and in the PbPb underlying event. The 
event selection for pythia+hydjet was based on reconstructed calorimeter jet information, as 
for the previous studies. 
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Figure 13: Distribution of the transverse momentum sum of tracks for three pj ranges, as a 
function of the distance AR to the leading and subleading jet axes. Results for the 0-30% cen- 
trality selection are shown for pythia+hydjet (upper row) and PbPb data (lower row). For 
each figure, the requirements on the dijet asymmetry Aj are given. Note that events with 
Aj > 0.22 are much rarer in the PYTHIA+HYDJET sample than in the data. Vertical bars are 
statistical and systematic uncertainties, combined in quadrature, the systematic contributions 
being 20%, independent of the bin. 



To derive the associated track spectrum for a given jet selection in data, the pj distribution 
of tracks inside a ring of radius AR = A(p'^ + Ar]^ and width of 0.08 around the jet axes 
was summed over all selected jets. The contribution of tracks from the underlying event, not 
associated with the jet, was estimated by summing the track pj distributions using an equal- 
size ring that was reflected around f] = 0, but at the same (p coordinate as the individual jet. For 
this procedure, jets in the region < 0.8 were excluded and only ring-radii up to AR = 0.8 
around the jet axes were considered, to avoid overlap between the signal jet region and the 
region used for background estimation. In addition, jets in the region \ fj\ > 1.6 were excluded 
to ensure the 0.8 radius rings would lie within the tracker acceptance. Statistical fluctuations in 
the underlying event limit this procedure to tracks with transverse momenta pj > 1 GeV/c. 

The summed px spectra from the jet regions and the underlying event regions were then sub- 
tracted, yielding the momentum distribution of charged tracks associated with the jets as a 
function of AR. 

The resulting distributions of associated track momentum as a function of track pj and AR are 



presented in Fig. 13 for four selections in dijet asymmetry, from Aj < 0.11 (left) to Aj > 0.33 
(right). For both data and pythia+hydjet results, the jet selections and Aj values are based on 
the reconstructed calorimeter jet momenta (Section |2.4| in order to have consistent event selec- 
tions for comparison. The middle bin boundary {Aj = 0.22) corresponds to the median of the 
Aj distribution for the 0-30% central PbPb events shown here. The top row shows the results 
for PYTHIA+HYDJET simulations. The track results shown for the PYTHIA+HYDJET simulations 
were found using the known ("truth") values of the track momenta from the embedded pythia 
events. The bottom row presents results for PbPb data. The track results shown for PbPb data 



3.3 Overall momentum balance of dijet events 
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were corrected for tracking efficiency and fake rates using corrections that were derived from 
PYTHIA+HYDJET simulations and from the reconstruction of single tracks embedded in data. 
In each panel, the area of each colored region in pj and AR corresponds to the total transverse 
momentum per event carried by tracks in this region. 

For the balanced-jet selection, Aj < 0.11, one sees qualitative agreement in the leading and 
subleading jet momentum distributions between pythia+hydjet (top) and data (bottom). In 
data and simulation, most of the leading and subleading jet momentum is carried by tracks 
with Pj > 8 GeV/c, with the data tracks having a slightly narrower AR distribution. A slightly 
larger fraction of the momentum for the subleading jets is carried by tracks at low pj and 
AR > 0.16 (i.e., beyond the second bin) in the data. 

Moving towards larger dijet imbalance, the major fraction of the leading jet momentum con- 
tinues to be carried by high-px tracks in data and simulation. For the Aj > 0.33 selection, it is 
important to recall that less than 10% of all PYTHIA dijet events fall in this category, and, as will 



be discussed in Section 3.3 those that do are overwhelmingly 3-jet events. 



While the overall change found in the leading jet shapes as a function of Aj is small, a strong 
modification of the track momentum composition of the subleading jets is seen, confirming the 
calorimeter determination of the dijet imbalance. The biggest difference between data and sim- 
ulation is found for tracks with pj < A GeV/c. For PYTHIA, the momentum in the subleading 
jet carried by these tracks is small and their radial distribution is nearly unchanged with Aj. 
However, for data, the relative contribution of low-px tracks grows with Aj, and an increasing 
fraction of those tracks is observed at large distances to the jet axis, extending out to AR = 0.8 
(the largest angular distance to the jet in this study). 

The major systematic uncertainties for the track-jet correlation measurement come from the 
^Independent uncertainty in the track reconstruction efficiency. The algorithmic track recon- 
struction efficiency, which averages 70% over the pj > 0.5 GeV/c and < 2.4 range included 
in this study, was determined from an independent PYTHIA+HYDJET sample, and from sim- 
ulated tracks embedded in data. Additional uncertainties are introduced by the underlying 
event subtraction procedure. The latter was studied by comparing the track-jet correlations 
seen in pure PYTHIA dijet events for generated particles with those seen in pythia+hydjet 
events after reconstruction and background subtraction. The size of the background subtrac- 
tion systematic uncertainty was further cross-checked in data by repeating the procedure for 
random ring-like regions in 0-30% central minimum bias events. In the end, an overall sys- 
tematic uncertainty of 20% per bin was assigned. This uncertainty is included in the combined 



statistical and systematic uncertainties shown in Fig. 13 



3.3 Overall momentum balance of dijet events 

The requirements of the background subtraction procedure limit the track-jet correlation study 
to tracks with pj > 1.0 GeV/c and AR < 0.8. Complementary information about the over- 
all momentum balance in the dijet events can be obtained using the projection of missing pj 
of reconstructed charged tracks onto the leading jet axis. For each event, this projection was 
calculated as 

A = 'Ej~Pt ( ~ ^Leading Jet ) , (2) 
i 

where the sum is over all tracks with pj > 0.5 GeV/c and \f]\ < 2.4. The results were then 
averaged over events to obtain (^f). No background subtraction was applied, which allows 
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Figure 14: Average missir\g transverse momentum, (^j), for tracks with pj > 0.5 GeV/c, pro- 
jected onto the leading jet axis (solid circles). The (^|) values are shown as a function of dijet 
asymmetry Aj for 30-100% centrality (left) and 0-30% centrality (right). For the solid circles, 
vertical bars and brackets represent the statistical and systematic uncertainties, respectively. 

Colored bands show the contribution to (^|) for five ranges of track pj. The top and bot- 
tom rows show results for pythia+hydjet and PbPb data, respectively. For the individual pj 
ranges, the statistical uncertainties are shown as vertical bars. 



this study to include the \rijet\ < 0.8 and 0.5 < pj^ < 1.0 GeV/c regions not accessible for the 
study in Section 3.2 The leading and subleading jets were again required to have \rj\ < 1.6. 



In Fig. 



14 



j) is shown as a function of Aj for two centrality bins, 30-100% (left) and 0-30% 



(right). Results for PYTHIA+HYDJET are presented in the top row, while the bottom row shows 
the results for PbPb data. Using tracks with < 2.4 and pj > 0.5 GeV/c, one sees that indeed 
the momentum balance of the events, shown as solid circles, is recovered within uncertainties. 



3.3 Overall momentum balance of dijet events 
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for both centrality ranges and even for events with large observed dijet asymmetry, in both 
data and simulation. This shows that the dijet momentum imbalance is not related to unde- 
tected activity in the event due to instrumental (e.g. gaps or inefficiencies in the calorimeter) or 
physics (e.g. neutrino production) effects. 
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Figure 15: Average missing transverse momentum, (^j), for tracks with pj > 0.5GeV/c, pro- 
jected onto the leading jet axis (solid circles). The values are shown as a function of dijet 
asymmetry Aj for 0-30% centrality, inside (AR < 0.8) one of the leading or subleading jet cones 
(left) and outside (AR > 0.8) the leading and subleading jet cones (right). For the solid circles, 
vertical bars and brackets represent the statistical and systematic uncertainties, respectively. 
For the individual pj ranges, the statistical uncertainties are shown as vertical bars. 

The figure also shows the contributions to (//|) for five transverse momentum ranges from 0.5- 
1 GeV/ c to Pt > 8 GeV/ c. The vertical bars for each range denote statistical uncertainties. For 

data and simulation, a large negative contribution to (|/|) (i.e., in the direction of the leading jet) 
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by the px > 8 GeV/c range is balanced by the combined contributions from the 0.5-8 GeV/c 
regions. Looking at the px < 8 GeV/c region in detail, important differences between data 
and simulation emerge. For pythia+hydjet both centrality ranges show a large balancing 
contribution from the intermediate pj region of 4-8 GeV/c, while the contribution from the 
two regions spanning 0.5-2 GeV/c is very small. In peripheral PbPb data, the contribution of 
0.5-2 GeV/c tracks relative to that from 4-8 GeV/c tracks is somewhat enhanced compared to 
the simulation. In central PbPb events, the relative contribution of low and intermediate-px 
tracks is actually the opposite of that seen in pythia+hydjet. In data, the 4-8 GeV/c region 
makes almost no contribution to the overall momentum balance, while a large fraction of the 
negative imbalance from high px is recovered in low-momentum tracks. 

The dominant systematic uncertainty for the px balance measurement comes from the px- 
dependent uncertainty in the track reconstruction efficiency and fake rate described in Sec- 
tion |3.2[ A 20% uncertainty was assigned to the final result, stemming from the residual dif- 
ference between the pythia generator-level and the reconstructed pythia+hydjet tracks at 
high Pj. This is combined with an absolute 3 GeV/c uncertainty that comes from the imperfect 
cancellation of the background tracks. The background effect was cross-checked in data from 



a random cone study in 0-30% central events similar to the study described in Section 3.2 The 



overall systematic uncertainty is shown as brackets in Figs.jMjand 15 

Further insight into the radial dependence of the momentum balance can be gained by studying 

(p'|.) separately for tracks inside cones of size AR = 0.8 around the leading and subleading jet 
axes, and for tracks outside of these cones. The results of this study for central events are 



shown in Fig. 15 for the in-cone balance and out-of-cone balance for MC and data. As the 
underlying PbPb event in both data and MC is not (^-symmetric on an event-by-event basis, 
the back-to-back requirement was tightened to A^i2 > 5n/6 for this study. 

One observes that for both data and MC an in-cone imbalance of (i/l) ~ —20 GeV/c is found for 
the A J > 0.33 selection. In both cases this is balanced by a corresponding out-of-cone imbalance 

of (p'|.) 20 GeV/c. However, in the PbPb data the out-of-cone contribution is carried almost 
entirely by tracks with 0.5 < px < 4 GeV/c whereas in MC more than 50% of the balance is 
carried by tracks with px > 4 GeV/c, with a negligible contribution from pj < 1 GeV/c. 

The pythia+hydjet results are indicative of semi-hard initial or final-state radiation as the 
underlying cause for large Aj events in the MC study. This has been confirmed by further 
studies which showed that in pythia the momentum balance in the transverse plane for events 
with large Aj can be restored if a third jet with px > 20 GeV/c, which is present in more than 
90% of these events, is included. This is in contrast to the results for large-Aj PbPb data, which 
show that a large part of the momentum balance is carried by soft particles (px < 2 GeV/c) and 
radiated at large angles to the jet axes (AR > 0.8). 



4 Summary 

The CMS detector has been used to study jet production in PbPb collisions at ^/s^ = 2.76 TeV. 
Jets were reconstructed using primarily the calorimeter information in a data sample corre- 
sponding to an integrated luminosity of Lint = 6.7 jih^^. Events having a leading jet with 
Px > 120 GeV/c and |//| < 2 were selected. As a function of centrality, dijet events with a 
subleading jet of px > 50 GeV/c and |;/ 1 < 2 were found to have an increasing momentum im- 
balance. Data were compared to PYTHIA dijet simulations for pp collisions at the same energy 
which were embedded into real heavy ion events. The momentum imbalances observed in the 
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data were significantly larger than those predicted by the simulations. While the relative im- 
balance between the leading and subleading jets increased with increasing collision centrality 
it was found to be largely independent of the leading jet pj, up to the highest pj region studied 
(f=i210 GeV/c). 

The angular distribution of jet fragmentation products has been explored by associating charged 
tracks with the dijets observed in the calorimeters. The calorimeter-based momentum imbal- 
ance is reflected in the associated track distributions, which show a softening and widening 
of the subleading jet fragmentation pattern for increasing dijet asymmetry, while the high-px 
components of the leading jet remain nearly imchanged. 

Studies of the missing transverse momentum projected on the jet axis have shown that the 
overall momentum balance can be recovered if tracks at low pj are included. In the PbPb data, 
but not in the simulations, a large fraction of the balancing momentum is carried by tracks 
having pj < 2 GeV/c. Comparing the momentum balance inside and outside of cones of 
AR = 0.8 around the leading and subleading jet axes demonstrates that a large contribution 
to the momentum balance in data arises from soft particles radiated at AR > 0.8 to the jets, a 
feature which is also not reproduced in PYTHIA calculations. 

In conclusion, a strong increase in the fraction of highly unbalanced jets has been seen in central 
PbPb collisions compared with peripheral collisions and model calculations, consistent with a 
high degree of jet quenching in the produced matter. A large fraction of the momentum bal- 
ance of these unbalanced jets is carried by low-px particles at large radial distance, in contrast 
to PYTHIA simulations embedded into heavy ion events. The results provide qualitative con- 
straints on the nature of the jet modification in PbPb collisions and quantitative input to models 
of the transport properties of the medium created in these collisions. 
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